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Abstract 

We compute the dimensionality dependence otn/s for charged black branes with 
Gauss-Bonnet correction. We find that both causality and stability constrain the 
\ value of Gauss-Bonnet coupling constant to be bounded by 1/4 in the infinite dimen- 

sionality limit. We further show that higher dimensionality stabilize the gravitational 
perturbation. The stabilization of the perturbation in higher dimensional space-time 
is a straightforward consequence of the Gauss-Bonnet coupling constant bound. 



: 1 Introduction 



The AdS/CFT correspondence [1-3] provides an interesting theoretical framework for 
studying relativistic hydrodynamics of strongly coupled gauge theories. The result of 
RHIC experiment on the viscosity /entropy ratio turns out to be in favor of the predic- 
tion of AdS/CFT [4-6]. Some attempt has been made to map the entire process of RHIC 
experiment in terms of gravity dual [7]. The way to include chemical potential in the 
theory was figured out in [8,9]. Phases of these theories were also discussed in [9-13]. 

It had been conjectured that the viscosity value of theories with gravity dual may give 
a lower bound for the rj/s = 7- for all possible liquid [14]. However, in the presence of 
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higher-derivative gravity corrections, the viscosity bound and causality are also violated as 
a consequence [15-22]. 

The higher derivative gravity terms are related to the (in)stability issues of black holes. 
The black hole stability issues are a crucial problem because black hole solutions are 
no longer unique in spacetime with higher than four dimensions. The instability of IP- 
dimensional asymptotically flat Einstein-Gauss-Bonnet black holes has been discussed by 
several authors [23,24]. Their results show that for the gravitational perturbations of 
Schwarzschild black holes in D > 5 Gauss-Bonnet gravity, the instability occurs only for 
D = 5 and D = 6 cases at large value of a' [24]. In the previous paper [20], we computed 
the charge dependence of rj/s for Gauss-Bonnet theory and noticed that charges introduced 
instability of the black brane even in the range < A < 0.09. 

The purpose of this paper is to perform a complete computation of rj/s including the 
charge and Gauss-Bonnet correction to an arbitrary dimensionality , and to determine the 
causality and stability constraints on the parameters of the black hole. Both the causality 
and stability constraints give the same result that A should be bounded by 1/4 for an 
arbitrary high dimensionality. We further find that higher dimensionality stabilize the 
tensor type perturbation. 

2 Viscosity to entropy density ratio 

We have explored the charge dependence of rj/s in the presence of Gauss-Bonnet term for 
five-dimensional AdS black branes [20]. In this section, we generalize the previous result 
on rj/s [20] to .D-dimensional cases. Let us start by introducing the following action in D 
dimensions which includes Gauss-Bonnet terms and U(l) gauge field: 

1 = 16^J dDx ^( R ' 2A+a ' ( R ^ R<XVPa - AR ^ V + R 2 )-^G D F^F^, (2.1) 
where a' is a (positive) Gauss-Bonnet coupling constant with dimension (length) 2 and the 
field strength is defied as F fll/ (x) = d fJt A l ,(x) — d l ,A fl (x). The thermodynamics and geometric 
properties of black objects in Gauss-Bonnet gravity were studied in several papers [25-31]. 
The charged black brane solution in D dimensions for this action is described by [28] 
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ds 2 = -H(r)N 2 dt 2 + H-\r)dr 2 + ^-darW, (2.2a) 
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where a and a' are connected by a relation a = (D — 4)(D — 3)a', A = a/l 2 , a = 
and the parameter I corresponds to AdS radius. The horizon is located at r = r H 
gravitational mass M and the charge Q are expressed as 
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Taking the limit a' — * 0, the solution corresponds to one for Reissner-Nordstrom-AdS 
(RN-AdS). The hydrodynamic analysis in this background has been done in [32,33]. 

The constant A^ 2 in the metric fl2.2al) can be fixed at the boundary whose geometry 
would reduce to flat Minkowski metric conformaly, i.e. ds 2 oc —c 2 dt 2 + dx 2 . On the 
boundary r — >• oo, we have 



H{r)N 2 ^ - 



I 2 



so that A^ is found to be 



N 1 



4A 



(2.4) 



Note that the boundary speed of light is specified to be unity c = 1. Eq. (12.41) implies that 
the significant value of A lies in the region A < 1/4. We will confirm this result from the 
causality and stability analysis in section 3 and 4. 
The temperature at the event horizon is defined as 

1 d^/gTt A^r^ 



T 



^y/grr dr 



((£> - 1) - (£> - 3)a) . 



(2.5) 



The black brane approaches extremal when a 
given by [27] 

1 



(i.e. T — > 0). The entropy density is 



~D-2 



We will calculate the shear viscosity of the boundary theory using the Kubo formula 

77 = lim ±- / dtdxe~ wt < [T xy (x),T xy (0)} >= - lim -lmG xy , xy (cu, 0), (2.7) 

where G(cu,0) is the retarded Greens function for T xy : 

G xy , xy (uj,k) = -ij dtdxe lk - M 6(t) < [T xy (x),T xy (0)\ 

It is convenient to introduce coordinate in the following computation 
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(2.9) 



We now study the tensor type perturbation h y (t,x 3 ,z) = (f)(t,x 3 ,z) on the black brane 
background of the form 

D-2 



ds 2 = -f(z)N 2 dt 2 + 



dz 2 



+ 



b 2 f(z) b 2 l 2 
where b = Using Fourier decomposition 

. r & D ~ l k 

<p{t,x 3 ,z) 



2<f>(t, x 3 , z)dxdy + ^ dx 2 J , 



i=i 



<j>{k,z), 



(2n) D - 1 

we can obtain the following linearized equation of motion for (f>(z) from the Einstein-Gauss- 
Bonnet-Maxwell field equation: 



g{z)<p" + g\z)<p' + g 2( f) = 



(2.10) 



where 



9(z) = z D - 2 f {1-^3 [z-V + z-\D - 5)/] } 
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and the prime denotes the derivative with respect to z. 



(2.11) 
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For the convenient calculation of the shear viscosity, we further introduce a new vari- 
able u{= -). Then we solve the equation of motion for transverse graviton Eq. fl2.10l) in 
hydrodynamic regime i.e. small to and k. The solution to Eq. fl2.10l) as 

</>(z) = (l-u) v F(u), (2.12) 

where F(u) is a regular function at the horizon u = 1, so that the singularity at the 
horizon might be extracted. The parameter v can be fixed as v = ±ia;/47rT by substituting 
Eq.f l2.12l) into the equation of motion. Usually we choose 

v = -i^— (2.13) 
4ttT v ; 

as the incoming wave condition. To obtain the shear viscosity via Kubo formula (12.71) . we 
only need know the u — > behavior of the transverse graviton, so it is sufficient to expand 
F[u) in terms of frequencies up to the linear order of u(= %^kTv\ 

F{u) = F (u) + vF x (u) + C(z/ 2 , k 2 ). (2.14) 

Expanding f)2.10p to the first order of we get the following form, 

[g(u)F'(u)}' - v (J—g( u ))'F(u) - -^g( u )F'(u) = 0. (2.15) 
\l — u J 1 — u 

Substituting the series expansion (I2.14p into the equation (12.151) . one can get the equations 
of motion for F (u) and Fi(u) recursively. Following the procedure given in [20], we easily 
get 

F (u) = C, (const.). (2.16) 

and 

= 7^7 + "7TV ( 2 - 17 ) 

l-u g(u) 

The integration constant C2 can be fixed by the regularity condition of Fi(u) at the horizon. 
So, the regularity condition at u — 1 implies 



2A 

{{D - 1) - (D - 3)a) (1 - p^((D - 1) - (D - 3)a)) C. (2.1S) 



C 2 = - 

The remaining constant C is estimated in terms of boundary value of the field 

lim <b{z 



so that we could fix 

C = <j>W(l + 0(u)\ (2.19) 

Now let us calculate the retarded Green function. Using the equation of motion, the 
action reduces to the surface terms. The relevant part is given as 

rf-'jV f d D - l k 



Near the boundary u = e, using the obtained perturbative solution for (f>(u), we can get 
((D - 1) - (D - 3)a) (1 - ^(p - 1) - (D - 3)a)) 



e = -u- 



Therefore we can read off the correlation function from the relation (12.81) , 

G xy xy (u, k) = -iujjr^r (jkz) (l " ^3 l(D -1)-(D- 3)o]) + 0(uA k 2 ), (2.22) 

where we subtracted contact terms. Then finally, we can obtain the shear viscosity by 
using Kubo formula (12.71) . 

" = (iS) ( l -w^v- D - 1) -^ D - 3)a ■ < 2 - 23) 

The ratio of the shear viscosity to the entropy density is found to be 

^s( 1 -z^3K B - 1 >- (B - 3 > a 0' (2 ' 24) 

The above result agrees with [20] when D = 5. Since A is bounded by 1/4, the shear 
viscosity never approaches zero in higher than 5D Gauss-Bonnet theory. When a = (no 
charges), rj/s = (1 — 2 ^ D 3 ^ )/ (47r), we recover the result in Ref. [16]. It is also worth 
noting that for extremal case (a = ;§E§), the ratio of the shear viscosity to entropy density 
receives no corrections from Gauss-Bonnet terms. In the next two sections, we will show 
explicitly causality and stability impose more rigorous constraint on the value of A and the 
upper bound of A is 1/4. 
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3 Causality constraints 



The authors in [16, 17] demonstrated that the causality could be violated if one introduced 
Gauss-Bonnet terms. In [20], it is found that the presence of charge does not contribute 
to causality and the result of [16, 17] is universal for charged black branes. In this section, 
we investigate the dimensionality dependence of the causality constraints. 

Due to higher derivative terms in the gravity action, the equation f 1 2 . X j) for the prop- 
agation of a transverse graviton differs from that of a minimally coupled massless scalar 
field propagating in the same background geometry. Writing the wave function as 

(f)(x 3 ,z) =e~ Mkz+ik3Xs , (3.1) 

and taking large momenta limit k^ — > oo, one can find that the equation of motion f 1 2 . X j) 
reduces to 

Wgf v ~ 0, (3.2) 

where the effective metric is given by 

,effj„)iJ™i' _ A7-2 /•/ \ / i.2 i 1 i 2 > ' 



dsi e = g^dx^dx" = N 2 f(r) ( -d? + ^da£J + — dr^. (3.3) 



Note that c 2 can be interpreted as the local speed of graviton: 

an _ N*f 1 ~ Xo=jo^ U" + (D-5)(D- 6KV + 2(D - 5)z^f) 

^ Z 2 l__2A_[ z -l / / + z -2 (D _ 5)/ 2 ] 



We can expand c 2 near the boundary z — > oo, 

(D 2 -5D + 10)(l + a) 



cl-1 



9 V. 2(D-3)(D-4) 

(D-l)(l + a) 1 + a 



(D - 3)(D - 4)(1 - 4A) 2v^4A7 z d ~ x + °^ ] ' ^ 
As the local speed of graviton should be smaller than 1 (the local speed of light of the 
boundary CFT), we require 

(D 2 -5D + 10) (D-l) 1 

' 2(£>-3)(£>-4) + (£>-3)(£>-4)(l-4A) ~ 2v^4A ~ °' (3 ' 6) 

The above formula leads to 

D 4 — 10D 3 + AID 2 — 92D + 96 

Acausality < 4(jD 2 _ 5D + 1Q )2 " ^.7) 
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without any charge dependence. As D is large enough, the above formula becomes 

lim A causa iit y < -■ (3.8) 

D^oo 4 

We can rewrite the above equation from the relation A = (D — 3)(D — 4)a'/l 2 , 

a' D 4 - 10D 3 + AID 2 - 92D + 96 

¥ - A(D 2 -5D + 10) 2 (D - 3){D - 4) ' ^ 
Figure 1 demonstrates that causality constrains the value of A. When D = 5, it goes as 
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Figure 1: The minimal value of A constrained by causality. The upper bound of A is 1/4. 

A < 0.09 same as the result of [17] and when D — > oo, A is bounded by 1/4. 
Now, we rewrite the wave function in a Schrddinger form, 

1 



( ' 2 " +^(z(r.))V>="V, P- 



dr 2 

where ip (z(r*)) and the potential is defined by 

z/> = K(z)<f>, K(z) 



dz Nf{z) 
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From the geodesic equation of motion 



9 » v ds ds ' 
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(3.10) 



(3.11) 



(3.12) 



and the Bohr-Sommerfield quantization condition 

1 

4 



dnJco 2 - k 2 c 2 g = (n- 7 )tt, (3.13) 



one can find that the group velocity of the test particle along the geodesic line is given 
by [17] 



v 9 = 7JT - c * > L ( 3 - 14 ) 



do; 

Therefore, signals in the boundary CFT propagate outside of the light cone and micro- 
causality violation happens (for a more detailed and explicit discussion on causality vio- 
lation, see [17]). Now we can conclude that as dimensions of space-time go up, causality 
restricts the value of A in the region A < 1/4. In next section, we will prove that in the 
extremal limit a — > -§r|, the stability of the black brane also requires that A should also 
be bounded by 1/4. 

4 Stability constraints 

In [20], it was found that apart from the causality violation, for RN-AdS black brane in 
Gauss-Bonnet theory, the charges give new instability of the black brane within the window 
of < A < 0.09. Now, we will show that higher D stabilize the gravitational perturbation. 

From Figure EJ we can see that the Schrodinger potential develops a negative gap near 
the horizon, but the gap becomes flatter as D increases. We will now show that while in 




Figure 2: Schrodinger potential V(u) as a function of u (u = -) and D for A = 0.24 and a = jy^- 

the large momentum limit, the negative- valued potential leads to instability of the black 

9 



brane, higher D tends to suppress those unstable perturbations. In the large momenta 
limit —>■ oo, the dominant contribution to the potential is given by k 2 c 2 . In [20], it was 
found that for near extremal cases, c 2 can be negative near the horizon and V ~ k 2 c 2 can 
be deep enough (see Figure [2]). Thus bound states can live in the negative- valued well. 
The negative energy bound state corresponds to modes of tachyonic mass on Minkowski 
slices [34] and signals an instability of the black brane [23,24]. As D goes up, we will see 
new physics in the following. Let us expand c 2 in series of (1 — -), 

c 2 =[(D - 3)a - D - 1] (l + v 7 ! - 4A) {D 2 [4A 2 (a - l) 2 + 2(a + 1)A - l] 
—D [8A 2 (3a 2 - 4a + 1) + 2A(a + 7)] + [A 2 (3a - l) 2 - 3A(a - 1) - 3] } 
{2{D - 4) [-3 + (2 - 6a)A + D(l + 2(a - l^)]}" 1 (1 - -) + 0((1 - -) 2 ). (4.1) 



z 



Since < a < jy^f, and < ~ < 1, c 2 will be negative, if 
{ J D 2 [4A 2 (a-l) 2 + 2(a+l)A-l] 

—D [8A 2 (3a 2 - 4a + 1) + 2A(a + 7)] + [A 2 (3a - l) 2 - 3A(a - 1) - 3] } 

{2(-4 + D) [-3 + (2 - 6a)A + D(l + 2(-l + a)A)]} _1 < 0. (4.2) 

From the above formula, we find the critical value of A, 

\ c (D,a)= ^{ - (D-l)(D-6) + (D-3)(D + 2)a 

+ j(£> - 1) 2 (5D 2 - 40D + 84) + (D - 3) 2 (5L> 2 - 24D + 52)a 2 

l _ 2 

-6a(D 2 -8D + 20)(D-3)(D-l)| 2 ||l + (D-3)a-D| . (4.3) 
Above the line of A c , c 2 can be negative (see figure [3]). The minimal value of A c can be 



obtained in the limit a 



■ D-l \ 



A ■ = j P-3)P-4) (44) 



That is to say 



a' 



JV Cimin i(D-l)(D-2Y (45) 
When D = 5, A c> m in = we recover the result obtained in ref. [20]. Eq (l4.4j) indicates that 
for any value of a, the quasinormal modes (QNMs) become stable under the line A Cj min . 
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Figure 3: The minimal value of A constrained by instability in the limit a — > 7533- The figure 
shows that A c is bounded by 0.25 

As the value of D increases, one finds that A c> m i n is also bounded by 1/4, i.e. 

lim A c = 7 (4.6) 

(£>,a)-(oc,2jEi) 4 

Note that this value is obtained in the extremal limit. Different from causality violation, 
the stability of the black brane depends on the charge. It would be very interesting to see 
for fixed value of charge, for which value of A the QNMs become stable. Eq (l4.3l) tells us 
that for A < X c (D,a), the black brane is always stable. Actually, as pointed out in [20], 
for fixed D the two lines A c (a) and A causa ii ty separates the physics into four regions in 
(a, A) space: consistent region; only causality violation region; only unstable modes region; 
causality violation and unstable modes region(see figure 4 in [20] for more details). In 
order to show explicitly the behavior of gravitational perturbation in higher dimensions, 
we solve the Schrodinger equation (13.101) with negative valued potential numerically and 
find some unstable QNMs (see tables OQ [2] and [3]) • Among these tables, we can find that 
the real part of to is vanishing, while the imaginary part of uj is positive. 

Table 1 demonstrates that the unstable modes of the black brane are suppressed as 
D increases. This confirms the result obtained in Ref. [24]. From table 2 and 3, we see 
that lower value of charge (a) and A stabilize the perturbation, while the lower value of 
D strengthens the instability. The reason for why higher D stabilize the perturbation is 
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n 


A — 20 


A — 18 


A — 1 6 


A — 1 2 


A — 10 


U 


1 99 797R? 


Q9 ^o^? 


fifi 001 9? 






7 




UZ .OO ± eft 








8 


49.7368i 


30.9208i 


16.8584 






9 


34.7918i 


20.0792i 


7.5717i 






10 


24.7602i 


11.9079i 









Table 1: Unstable QNMs for charged GB black brane perturbation of tensor type for fixed charge 
(a = 1.20) and k 3 = 500. 



D 


a = 1.4 


a = 1.2 


a = 1.0 


a = 0.8 


a = 0.6 


a = 0.4 


6 


129.0012 


122.7276i 


110.9159^ 


91.1455i 


59.3517i 


7.0564i 


7 


79.3218i 


75.3435? 


66.5505i 


37.3177i 


24.3269i 




8 


70.0214i 


49.7368i 


42.0717i 


26.6617i 






9 




32.6603i 


27.1475i 


12.6057i 






10 




24.7602i 


20.6077i 


4.6077i 







Table 2: Unstable QNMs for charged GB black brane perturbation of tensor type for fixed A 
(A = 0.20) and k$ = 500. Note that a = 1.4 exceeds the maximal value of charge permitted for 
9- and 10-dimensional charged black brane and thus we leave the frequency blank there. 



A 


a = 1.6 


a = 1.4 


a = 1.2 


a = 1.0 


a = 0.8 


a = 0.6 


a = 0.4 


0.20 


131.38692 


129.0012 


122.7276i 


110.91592 


91.14552 


59.3517i 


7.0564i 


0.15 


66.3475i 


62.7713i 


53.5830i 


36.9280i 


10.0682i 






0.10 


16.0152 


10.3104i 












0.08 


0.9712i 















Table 3: Unstable QNMs for charged GB black brane perturbation of tensor type for fixed 
dimensionality (D = 6) and ks = 500. 
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because that no matter how big D is, A (i.e. A — (D — 3)(D — 4) a' /I 2 ) is bounded by 1/4 
which means that for fixed I, a' — > as D increases. Moreover, a' — > corresponds to 
vanishing Gauss-Bonnet correction in ( 12. ip and charged black branes in that regime are 
stable. 

5 Conclusions and discussions 

In summary, we have computed the dimensionality dependence of rj/s for charged black 
branes with Gauss-Bonnet correction. The ratio of the shear viscosity to entropy density 
in D-dimensional space-time was found to be r]/s = 4- (l — -£^[(D — 1) — (D — 3)a]). 
When D = 5, we can recover the result found in [20]. It is worth noticing that for non-zero 
charge the viscosity can never approaches even in 5D case. 

While in [35], it was always assumed that A < 7, in this paper we have shown explicitly 
that both causality and stability constrained the value of A to be bound by 1/4 in the limit 
D — > 00, but for fixed D, these two constraints are different. It is interesting to notice that 
from different physical processes (causality and stability), we obtain the same bound. One 
may further check whether this is a coincidence or not. 

The instability of charged black brane with Gauss-Bonnet correction was also analyzed 
in this paper. The result shows while higher value of charge (a) and A strengthen the 
perturbation, the unstable modes of charged black brane are suppressed as D increases. 
The suppression of the unstable modes for higher D can be explained from the fact that A 
is bounded by 1/4 and a ' approaches zero as D increases. 
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